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Abstract. 

The ability to monitor nutrient and other environmental conditions with high 
sensitivity is crucial for cell growth and survival. Sensory adaptation allows a cell 
to recover its sensitivity after a transient response to a shift in the strength of 
extracellular stimulus. The working principles of adaptation have been established 
previously based on rate equations which do not consider fluctuations in a thermal 
environment. Recently, G. Lan et al. (Nature Phys., 8:422-8, 2012) performed a 
detailed analysis of a stochastic model for the E. coli sensory network. They showed 
that accurate adaptation is possible only when the system operates in a nonequilibrium 
steady-state (NESS). They further proposed an energy-speed-accuracy (ESA) trade¬ 
off relation. We present here analytic results on the NESS of the model through a 
mapping to a one-dimensional birth-death process. An exact expression for the entropy 
production rate is also derived. Based on these results, we are able to discuss the ESA 
relation in a more general setting. Our study suggests that the adaptation error can 
be reduced exponentially as the methylation range increases. Finally, we show that a 
nonequilibrium phase transition exists in the infinite methylation range limit, despite 
the fact that the model contains only two discrete variables. 


1. Introduction 

As a paradigmatic example of environmental monitoring in biology, the E. coli 
chemotactic sensory system has been studied extensively over the years p]. Its core 
component is the transmembrane methyl-accepting chemotaxis protein (MCP) receptor. 
MCP binds selectively to ligands outside the cytoplasmic membrane and modulates 
the activity of its downstream signal transduction pathway in a way that depends on 
its methylation state. Two aspects are recognized to be crucial to the performance 
of the sensory network in the biological context: sensitivity of detection in a noisy 
environment, and adaptation to maintain that sensitivity over a broad range of ligand 
concentrations. With regard to high sensitivity to diffusing chemicals in the surrounding 
medium, Berg and Purcell |2J presented an optimal strategy in 1977 based on simple 
physical considerations. They showed that the measured chemotactic sensitivity of 
E. coli approaches that of the optimal design. Further indication of the organism’s 
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optimal performance is found in its nearly perfect adaptation over five decades in ligand 
concentration. The latter property is shown to hold even when proteins on the sensory 
network are expressed away from their natural levels [3j. To explain this remarkable 
behavior, Barkai and Leibler (BL) [3] introduced a simple model where the MCP 
methylation/demethylation rates are linked to the downstream activity. The system 
reaches a steady state only when its activity is at the level required by the balance of 
methylation and demethylation currents, ft was soon pointed out by Yi et al.|5] that 
the BL scheme is in effect implementing an integral feedback control which is widely 
used in engineering systems to achieve robust adaptation. Furthermore, an exhaustive 
search by Ma et al. [5] to identify all possible 3-node adaptive networks found integral 
feedback control as one of only two core motifs that enable perfect adaptation. 

In a separate development, there have been much progress in recent years in 
understanding fluctuation phenomena in non-equilibrium systems whose dynamics do 
not satisfy detailed balance HUH El m ED H2]. One particular aspect of fluctuations 
in a nonequilibrium steady state (NESS) is the production of system’s entropy and its 
subsequent release as heat to the environment [13]. In this respect, the generic behavior 
of nonequilibrium systems studied in the statistical physics community is shared by 
molecular processes in a living cell. A well-known example is the kinetic proof-reading 
discussed by J. J. Hopfield in 1974 [T4]. Here, the molecular machinery to carry out 
DNA replication can achieve a much lower error rate by operating out of equilibrium. 
It may be argued that employing energy flux to enable or improve the performance of 
a molecular circuit is a common practice in biology. However, there have been only a 
few examples so far where the details are convincingly elucidated da m nzi ns]. 

The integral feedback control for adaptation requires asymmetric interactions 
between the output node and the integration node, which can only be realized by 
systems in a NESS. The issue of energy cost to maintain such a state was addressed in 
a recent study by Lan et al. [19]. One of their main findings is a relation among the 
energy dissipation rate, adaptation speed and adaptation accuracy (ESA), which they 
suggested to hold generally. Their result is based on a sensory network model which 
has been shown to reproduce most of the experimental data on the MCP receptor in E. 
coli dl- 

Despite its intuitive appeal, the ESA relation has not been derived from the more 
general results in the literature regarding the NESS. Should there be a fundamental 
connection between the adaptation accuracy and energy dissipation rate? Specifically, 
is there a lower bound for energy dissipation rate to achieve a given adaptation accuracy? 
To clarify this and other issues, it is necessary to perform a more comprehensive study 
of the sensory network model. Due to the conceptual importance of the sensory network 
model, a rigorous discussion is desirable. 

The paper is organized as follows. The biological background of adaptation and 
the model by Lan et al. are introduced in Section [2j followed by a detailed analysis of 
the NESS of the system in Section [3] In Section [4] we drive an exact expression for the 
energy dissipation rate and compare it with the ESA tradeoff relation. A nonequilibrium 
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phase transition of the system in the infinite methylation range limit is identified in 
Section [5] and its properties discussed. Section [6] contains a summary of our results and 
conclusions. Mathematical details of an approximate treatment of the NESS distribution 
is relegated to Appendix A. 

2. A model for sensory adaptation 

Here, we briefly introduce the transmembrane methyl-accepting chemotaxis protein 
(MCP) receptor, which is the core component for receiving signal and exercising 
adaptation [I]. MCPs regulate the clockwise-counterclockwise rotational switch of 
downstream flagellar motors which drive the run-and-tumble motion of an E. coli cell. A 
simple cartoon of this receptor is illustrated in Figure [l]( a). The activity of the receptor 
can be described by a binary variable a: a = 1 for the active state and a = 0 for 
the inactive state. The transition rate between the two states depends on the external 
ligand concentration (i.e., signal strength) s and the internal methylation level m. m 
ranges from 0 to m o, with mo = 4 for a single MCP. The methylation level can be 
increased by enzyme CheR and decreased by enzyme CheB, in a way that depends on 
the activity of MCP. Figure [ljb) illustrates response of the MCP receptor to a stepwise 
signal s obtained from experimental measurements. The mean activity (a(f)) changes 
sharply in a short time window r a less than a second, and recovers slowly over a much 
longer time scale r m , of the order of a minute, due to the slow change of average internal 
methylation level ( m(t )). The output recovery after a transient response to external 
stimuli is called adaptation. The performance of adaptation is characterized by the 
adaptation error which can be defined as the ratio between the final shift of activity and 
the relative change of signal strength, as illustrated in Figure [l](b) . 




Figure 1 . (a) The MCP receptor, (b) Mean response of the MCP receptor to a 
stepwise signal. 

A Markov network model with internal states specified by (a, m ) was proposed for 
a single MCP by Lan et al. [T9] , as shown in Figure [2] Transition between active and 
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inactive states at a given methylation level takes place on the time scale r a , with rates 
given by 

uii(m,s) = —expf—A E(m,s)\ u ; 0 (m, s) = —expf— — AE(m, s)Y 

r a \ 2 ) r a \ 2 / 

where (3 — \/T is the inverse temperature, and 
A E(m, s ) = eo(mi — m) + /(s) 

the free energy difference between states (0, m ) and (1, m), with /(s) = ln[(l+s//ij)/(l+ 
s/K a )\. Here e 0 > 0 is the methylation energy, m\ an offset methylation level [2U1 IT], 
and Ki and K a (^> K L ) are equilibrium constants for ligand binding to the receptor in the 
inactive and active states, respectively. Transition between different methylation levels 
takes place on the time scale r m , with rates indicated in Figure [2] when the receptor is 
inactive, the rate of methylation (assisted by enzyme CheR) is Ken while the rate of 
demethylation is olK C r ; when the receptor is active, the rate of demethylation (assisted 
by CheB) is Kcb while the rate of methylation is aKcB■ An estimate of the typical 
methylation/demethylation cycle time is given by T m = k^ B + k^ l R . 

The parameter a specifies the degree of disequilibrium in the system. When 
ot = cxeq = exp(/3e 0 /2), global detailed balance condition is satisfied, in which case the 
Boltzmann distribution governed by a free energy function is recovered. For a < oleq-, 
the system is driven out of equilibrium with generally different properties which we 
study using both analytical and numerical methods. Therefore a describes the strength 
of driving to keep the receptor to operate under out of equilibrium conditions. 

In the numerical examples presented below, we adopt the parameter values as 
suggested in Ref. [T9] : m\ = 1, K t = 18.2/iM, K a = 3000/uM, (3—1 (kT as 
the unit of energy) and eo = 2. The time constants are chosen as r a = 0.1s, and 
kcB = kcR = 0.01s -1 . For this parameter set, = e. To simplify the notation, we 
write uj\(m, s),uo(m, s) as oj\(m) and cuo(m) respectively. 
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Figure 2. The Markov network model of a single receptor MCP in E. coli. Red arrows 
indicate existence of a futile cycle at small a, which is essential for adaptation. 
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3.1. Condition for adaptation 


We first revisit the condition for adaptation first obtained by Lan et al. [19]. For the 
model introduced in Sec. 2, the master equation for the joint probability P(a, m) takes 
the form, 


dP( 0, m) 
dt 

dP( 1, m) 
dt 


= k C RP{ 0, m — 1) + ak C RP( 0, m + 1) + cui(m)P(l, m) 

- [k C R + ak CR + u o (m)]P(0,m), 

= a/cc- B P(l, m — 1) + k CB P( 1, m + 1) + ca 0 (m)P(0, m) 

- [akcB + k C B + uJi(m)\P(l,m). 

From the above, we obtain the evolution equations for the moments 
E a ,m mP ( a > m ) and («) = 
d(a) 


dt 

d(m) 


^uj o (m)P(0,m) — LOi(m)P(l, m) 


— (1 — a)(k C R + k C B) (~{a) + v - ~r —) 

V kru + kn r/ 


+ B L 


(la) 

(lb) 
m) = 

(2 a) 

(26) 


dt V kcR + kcB' 

Here Pi = akcnP( 0,0) + kcBP( 1,0) — kcRP(0,mo ) — afccBP(l, mo) depends on the 
probabilities for the extreme methylation states m = 0 and m = m 0 . 

In a steady environment of constant ligand concentration s, the system is expected 
to reach a steady state in a time r m where both (m) and (a) assume constant values. 


Setting the right-hand-side of Eq. (2b) to zero yields, 

kcR 


(a) a s 


+ 


P. 


(3) 


+ kcB 1 — a kcR + /cob 
S ince the methylation and demethylation rates kcR and kcB are assumed to be constants 
in the model, the first term a 0 = kcR/(kcR + k C B ) on the right-hand-side of Eq. ([3]) is 
independent of s. Figure [3] shows a s against s for three different values of a, obtained 
from numerically exact solution of the model in the NESS. The steady-state activity a s 
is centered around a 0 (dashed line) over a large range of s for a < 1, but not so for 
a > 1. The “adaptation error” 

1 Pi 

6 — \CL S CLq\ — 


(4) 

1 — ol kcR + kcB 

is essentially controlled by the size of the boundary term Pi. For a < 1, Pi is small over 
a broad range of s. As we shall see in the next section, the NESS distribution in this 
case is indeed centered in the middle of the allowed methylation range. This is however 
not the case when a > 1. 

The transient response to a signal ramp also exhibits qualitatively different behavior 
for a < 1 and a > 1. Figure [4] shows our results obtained by numerically integrating 


the master equations (la) and (lb) at three different values of a, upon a jump in ligand 


concentration from 10 K t to 15A,; at t — 0. The initial response of (a) to the signal 
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Figure 3. Mean receptor activity against ligand concentration for three different 
values of the nonequilibrium parameter a. The dash line indicates the value a o- Here 
the methylation range mo = 4. 


ramp is qualitatively similar in the three cases, i.e., a fast depression of receptor activity 
to a near plateau value in a time of order r a . However, opposite behavior is seen at 
longer times, in concert with the change in methylation level as seen in Fig. |4])b). For 
a > 1, a further decrease of the mean activity is seen when the methylation level starts 
to decrease in response to the change in (a). On the other hand, when a < 1, the 
methylation level increases in accordance with Eq. (2b), eventually restoring the mean 
activity to a value close to the pre-stimulus level. The latter is precisely the scenario for 
adaptation that employs a change in the methylation level to offset the activity change 
effected by the shift in signal strength. 

In summary, both the steady-state activity and the transient response to a signal 
ramp show qualitatively different behavior below and above a = 1. We thus conclude 
that the condition for adaptation in this model is a < 1. 


3.2. The NESS distribution 


In the previous subsection, we obtained the condition for adaptation by considering the 
moment equations with the help of numerical integration of the master equation. To 
gain a complete understanding of the NESS, it is necessary to calculate the distribution 
function P(a,m). Fortunately, for the model in question, this can be done under the 
“fast equilibrium” approximation facilitated by the separation of the time scales r a and 
r m r a . Then, P(l,m) and P(0,m) satisfy the local detailed balance 


P(l, m) 
P( 0, m) 


uJojm) 


+ 0(r a /r m ), 


(0 < m < mo). 


(5) 


Let P(m) = P(0,m) + P(l,m), we obtain, 

1 


P(0, m) = 


1 + exp[— f3AE(m, s)] 


P(m), P(l,m) — HPl ~f» l !’ j)l / W.(6) 

v h v ’ ; 1 + exp[-(3AE(m,s)] y J J 
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(a) 


(b) 


Figure 4. An initial steady state at s = lOA'* is perturbed by shifting the ligand 
concentration to s = 15A^ at t = 0. Results at three different values of a are shown: 
(a) change in the mean receptor activity ( a/t )) against f; (b) change in the methylation 
level (m(t)) against t. The methylation range m 0 is set to be four. 


With the help of ([5]), Eqs. ( JTa[ ) and (16) combine to yield 
dP(m) 


dt 


= b/m — 1 )P/m — 1) + d/m + l)P(m +1) — \b/m) + d/m)}P/m).( 7) 


Equation ([T]) defines a one-dimensional birth-death process with the birth and death 
rates given respectively by, 


b/m) = 


k CR + ak C B exp[-/3AE(s, m)] 


1 + exp[— {3AE(s, m)] 

Its steady-state distribution takes the form, 

b/m 


d(m) = 


ak C R + k C B exp[-/3AE(s, m)} 
1 + exp[— (3AE(s, m)] 


P(m + 1) = 


d/m + 1 ) 


P/m ) = P(0) Yl 


b(i) 


i =0 


d(i + 1)' 


( 8 ) 


Together with Eq. (|6j) the full NESS distribution is obtained. 

Consider the range of ligand concentrations where the receptor is functional, i.e., 
exp[— AE(m, s)] C 1 at m = 0 (inactive state favored) and exp[— AE(m, s)] 1 

at m = m 0 (active state favored). Consequently, the ratio b(m)/d(m + 1) changes 
monotonically between the limiting values 1/a and a as m increases from 0 to m q. 
Let m* be the value of m where b(m*)/d(m* + 1) ~ 1. According to Eq. ([8]), this 
is the methylation level where P(m) varies slowest with m, i.e., the stationary point 
of the distribution. For a < 1, b(m)/d(m + 1) > 1 on the low methylation side 
(m < m*) while b(m)/d(m + 1) < 1 on the high methylation side (m > m*). Therefore 
P/m) reaches its peak value at m*. The opposite situation happens for a > 1, where 
P/m ) initially decreases with m on the low methylation side, reaches its minimum 
value at rn* , and increases on the high methylation side. At a = 1, b/m) = d/m) 
so that P/m) = P(0)6(0)/d(m) becomes essentially flat especially when k C R = kcB- 
The general behavior of the NESS distributions in the two regimes are illustrated in 
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Figure [5j A more complete discussion of the functional form of these distributions and 
their dependence on s at different a values can be found in the Appendix A. 
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(a) (b) 

Figure 5. The steady-state distributions P(l,m), P(0,?n) and P{m) at too = 4, 
s = 10 Ki. (a) a = 0.1; (b) a = cyeq > 1. 


In the case a < 1, the signal level affects the shape of the distribution by shifting 
its peak position m* = m 1 + /(s)/e 0 — (/^eo)” 1 ln(fc CB //c ( y i? ) (see Appendix A) which 
coincides with the mean methylation level. This is a general feature for adaptation 
achieved through integral feedback control, i.e., the effect of external signal change is 
absorbed by a shift in the average methylation level. 


3.3. Adaptation error 


According to Eq. ([3]), the mean receptor activity (a) in the NESS depends on the signal 
level s only through the probabilities for the extreme methylation states at rri — 0 
and m = mo- For a < 1, P(m ) decreases rapidly away from the peak position at 
m*. Except very close to a = 1 which requires a separate treatment, the adaptation 
error as defined by Eq. (Jd]) can be estimated by the largest term in the expression for 
Bi. Let rri c j = min{m*,mo — m*} be the distance between m* and the closest extreme 
methylation state. With the help of Eq. (A.3), we obtain, 

— —m 2 d , rn d < ^^lna; 

In 2 a, otherwise. 


In e ~ 


l+a 2 "dt 

m d In a + ^ 


(9) 


2/3e 0 


Equation ([9]) shows that the adaptation error can be decreased by either increasing 
the methylation range m 0 or decreasing the parameter a that brings the system further 
away from equilibrium. At a given a < 1, increasing mo allows a greater functional 
range of the receptor and consequently larger values for m d , resulting in an exponential 
decrease of e. On the other hand, at a fixed mo, decreasing a increases the rate of 
exponential decay of e. However, when a is below a m = exp(—^eom 0 /2), the error 
basically saturates to a value bounded from below by e m ~ exp(— /3e 0 ml/8). These 
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observations are in agreement with the trends see in Fig. |6^a) where e is plotted against 
a for several different values of 


4. Energy dissipation and the ESA trade-off 

The nonequilibrium methylation/demethylation dynamics of the MCP receptor requires 
energy input [19]. Within the adaptation model considered here, the rate of energy 
dissipation can be calculated using the standard formula 


w = X'£j(x\x')i a u{xlx ' ) 


2/3 


X,X' 


u{X'\X)' 


( 10 ) 


Here oj(X\X') is the transition rate from state X' to state X, J( X\X') = 

u:(X\X')P( y X')—u:(X'\X)P( K X) is the net flux from X' to A", and P(X) is the probability 
for state A". For our purpose, it is convenient to rewrite the above equation in terms 
contributions from directed “elementary cycles” m An elementary cycle is a loop 
formed by nodes and edges of the network that cannot be further decomposed into 


smaller loops. Denoting by Q the Zth elementary cycle on the network, Eq. (10) can be 
rewritten as 


w = 1 £ J(C,)A(C,), 


( 11 ) 

and A(Ci ) = 


where J(Ci) is the probability flux associated with cycle Ci , 

£ eeCi \yl[oj{X'\X) / (jj(X\X')\, summed along the cycle. 

For the network model shown in Fig. [2j we define the mth elementary cycle to 
be the rectangle between methylation levels m and m + 1, directed counter-clockwise 
as indicated by the red arrows. It is simple to verify that the thermodynamic force 
A(Ci ) = 2\i\[q>eq/ a] is the same for all cycles. The cycle flux J{m ) = kcBP(l,m + l) — 
ak CB P(l, rn) can also be read off easily from the figure. From Eq. 0 we then obtain, 

1710 — 1 


2k 


In ——) ^ [fc CB P(l, m + 1) - ak C BP( 1, rn)] 


a 


CB 


In 


m =0 
Ot-EQ 


[(1 - a)a s - P( 1, 0) + aP(l, m 0 )\. 


(3 \ a 

With the help of Eq. (|3j) , we obtain finally the following exact expression for the energy 
dissipation in the NESS, 

2 


W = 


/3r n 


1 — a + I?2 ) In 


Ot-EQ 

a 


( 12 ) 


where B 2 = aP( 0, 0) — P(l, 0) + aP(l, m 0 ) — P(0, m 0 ) is a boundary term. 

Figure |6](b) shows W against a for selected values of m^. In all cases presented, a 
logarithmic increase on the far-from-equilibrium side (i.e., a -C 1) is seen, in agreement 
with Eq. (12). Dependence of W on m^, which enters only through the boundary term 
B 2 , is essentially negligible. This behavior can be understood from the fact that most 
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(a) (b) 

Figure 6. (a) Adaptation error and (b) energy dissipation over the methylation 

time scale r m against In a. Here the ligand concentration is chosen such that 
irid = m* = TOq/ 2 — 1 with =4,8 and 12. 


of the dissipation takes place in the loop centered around the peak position m* of the 
NESS distribution P(m). 

In Ref. [19], based on approximate solutions of the adaptation model at m 0 = 4, 
Lan et al. proposed the Energy-Speed-Accuracy tradeoff relation, 

W ~ (coal)u m ln(e 0 /e) (13) 


to capture the increase in energy dissipation to achieve higher accuracy of adaptation 
as a is reduced. Here cocx^ sets the appropriate energy scale for the problem, and 
(jj m = r m l . Comparing with our results Eqs. ([9]) and (12), we see that Eq. (13) needs to 
be modified to take into account the dependence of e on the distance ro^ from the actual 
mean methylation level to the boundaries of the full methylation range, i.e., m — 0 and 
m = m 0 . In addition, as we see from Fig. [6] , the adaptation error saturates to a value 
of the order of e m set by m o when a falls below a m , while the energy dissipation rate W 
keeps increasing. From the calculations presented above, we see that e is controlled by 
the probabilities for the rare events where the extreme methylation states are visited, 
while W is not sensitive to the actual methylation level itself. 


5. Phase transition 

As we have seen in Sec. 3, there is a qualitative change in the shape of the NESS 
distribution P{m ) at a = 1. For a > 1, P{m) is bimodal with peaks at the two ends 
of the methylation range from 0 to mo- The relative weight of the peaks is controlled 
by the signal strength s. On the other hand, for a < 1, P{m ) has a single peak 
in the middle of the methylation range. As the signal strength s varies, the peak 
position shifts accordingly but its shape remains more or less the same until either end 
of the methylation range is reached. As discussed previously by Lan et ah [19], the 
latter feature is crucial for the implementation of precise adaptation. In this section we 
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examine the transition between the two regimes in further detail with the help of the 
exact solution for the NESS. 

Let us first examine the behavior of the energy dissipation rate in the NESS given 
by Eq. (12) in the limit ro^ = min{m*,mo — m*} —> oo. For a < 1, Eq. (A.3) 


shows that the boundary probabilities vanish in this limit, and hence B 2 = 0. On the 
other hand, for a > 1, this limit implies the activation energies AE(0,s) —> +oo and 
AE^mo, s ) —)■ — oo. Therefore the receptor is nearly exclusively in the inactive state 
when the methylation level is close to zero, and exclusively in the active state when the 
methylation level is close to full. Then, the elementary loop current J{m) = 0 for all 
m which in turn yields vanishing dissipation. Here we encounter an interesting example 
where the detailed balance is violated by the kinetic rates but no dissipation actually 
takes place due to vanishing loop currents. Summarizing, we have in the limit —>■ oo, 


Wn = 


^(l-a)ln(oW«), 

0, 


0 < a < 1 ; 

1 < a < oieq- 


(14) 


The singular behavior of W against a indicates a true nonequilibrium transition in the 
model where the methylation range is infinite. 

According to Eq. §■ the adaptation error e can be made arbitrarily small in the 
entire adaptive phase a < a c = 1 by increasing m^. On the other hand, W can be 
made arbitrarily small at the same time by choosing an a close to a c . The energy 
dissipation is necessary to generate adaptive behavior, however, there does not appear 
to be a minimal value for the dissipation rate to support an arbitrarily accurate adaptive 
system. 

For a system with a finite methylation range, transition between the two phases 


is more gradual than what is described above. From Eqs. (A.3) and (A.4), one may 
identify a “correlation length” A ~ l/|lna| ~ jl — a\~ l . For m ( i > A, Eq. (14) can 


be directly applied. Corrections need to be considered when < A, based on exact 
results derived in previous sections. 


6. Conclusions 

In this paper, we report a detailed analytical study of the stochastic network model 
proposed by Lan et al. shown previously to describe well sensory adaptation in E. 
coli. To understand this system, we first derive moment equations which are closely 
related to the rate equations traditionally used to model this type of biological processes. 
The moment equations implement an integral feedback control scheme at the heart of 
the adaptive behavior. Adaptation in the model is achieved when the nonequilibrium 
parameter a < a c , where a c = 1 is less than its value cxeq when detailed balance is 
observed. By mapping the original ladder network to a one-dimensional birth-death 
process under the assumption of timescale separation, we obtain analytic expressions 
for the NESS distribution with qualitatively different behavior for a > 1 (nonadaptive 
phase) and a < 1 (adaptive phase). With the help of the exact results on the NESS 
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distribution, we compute the mean receptor activity from which its dependence on the 
external signal strength is obtained. In the adaptive phase, the adaptation error, which 
measures the deviation of the mean receptor activity from a suitable reference value, is 
found to decrease when the system is driven further out of equilibrium by reducing a, 
but approaches a saturated value for a < a m . Interestingly, at a given a , the adaptation 
error decreases exponentially with the number of the available methylation states before 
the extreme methylation levels are reached. 

We also derive an exact formula for the energy dissipation rate by using a cycle- 
decomposition technique. The energy dissipation rate is found to be insensitive 
to the size of the methylation range and also to timescale separation. Although 
our results confirm qualitatively the statement that adaptation within the molecular 
construct that implements integral feedback control requires nonequilibrium driving, 
there does not appear to be a lower bound on energy dissipation to achieve a given 
level of adaptation accuracy, in contrary to the Energy-Speed-Accuracy tradeoff relation 
proposed previously by Lan et al. 

Although the methylation range of a single MCP receptor is four, we have 
investigated the behavior of the system with arbitrary methylation range, especially 
when the methylation range mo is large. The extension allows us to examine various 
theoretical issues quantitatively. In the limit mo —)• oo, a true nonequilibrium phase 
transition at a = 1 can be identified. 

The current analysis only focuses on the static properties of the system. However, 
the transient response at short times is also an important component of the molecular 
adaptive circuit. An understanding of the adaptive behavior in a general setting based 
on thermodynamic principles is still lacking. In this respect, lessons may be drawn from 
recent developments in information thermodynamics [18, .25, [26] by viewing the adaptive 
circuit as an information processing machine. 
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Appendix A. Approximate expression of the NESS distribution 


In this Appendix we derive an approximate analytic expression for the NESS distribution 
given by Eq. ([8]). Taking the logarithm of the equation, we obtain, 


P(m + 1) 

In-—— 

P(0) 


m 


S ln 


b (i) 

d(i + 1) 


In 


d(0) 

d{m + 1) 


i=0 


(A.l) 


where 


4>(x) = In 


1 + ae x 


a + e x 
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and x(m) = —(3AE(s,m ) + ha.(kcB/kcR) — ( m — m*), with m* = mi + f(s)/e o — 
(,5eo) _1 \n{kcB/kcR)- It is straightforward to verify that cj)(x) is an odd function of x. 
The function (j>(x) is well approximated by a piece-wise linear function 


i]){x) 


— In a, 

azl T 

a+l X ’ 

In a, 


x < -£; 
-£ < x < £; 
x > £. 


(A.2) 


which has the same slope at x = 0 and same asymptotic values as x —> ±oo. Continuity 
requires the choice £ = In a. The two functions match each other well except near 
x = ±£. 

For a < 1, the peak of P{m ) is centered at m*. It is thus convenient to use P(m*) 
as the reference. Approximating the sum in Eq. (A.l) by an integral over ip{x), we 
obtain, 


In 


P(m ) 
P{m*) 




m*) In a + In -E- i n 2 ( 

' a(m) 1 —a 2peo 

« j/_* \ 


(m 


. m*P 


m < m_; 
m_ < m < m + ; 
m > m + . 


(A.3) 


Here m± = m* ± £/(/3eo). Since d{m ) has only a relatively weak dependence on m, we 
see that P(m ) is essentially gaussian within the interval (m_,m+), but turns to simple 
exponential decay outside the interval. 

For a > 1, P(m) achieves its minimum value at m*. In the neighborhood of the 
methylation boundaries, we have 


P{m ) ~ 


a m P(0)d(0)/d(m), m < m_; 

a m - m °P(mo)d(mo)/d(m), m > m + . 


(A.4) 


Noting that cf>(x) is an odd function of x, we have approximately P(m 0 ) — 

a m °- 2m *P(0)d(0)/d(m 0 ). 
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